Weakly interacting two-dimensional system of dipoles: limitations of mean-field theory 
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We consider a homogeneous 2D Bose gas with repulsive dipole-dipole interactions. The ground- 
state equation of state, calculated using the Diffusion Monte Carlo method, shows quantitative 
differences with predictions of commonly used Gross-Pitaevskii mean-field theory. The static struc- 
ture factor, pair distribution function and condensate fraction are calculated in a wide range of 
the gas parameter. Differences with mean-field theory are reflected in the frequency of the lowest 
"breathing" mode for harmonically trapped systems. 
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The study of quasi-two-dimensional Bose gases at 
ultra-low temperatures has become a very active area 
of research. The role of correlations and quantum fluc- 
tuations is greatly enhanced in reduced dimensional- 
ity making a two-dimensional (2D) system well suited 
for studying beyond mean-field effects. The superfluid- 
normal phase transition occurs at a finite-temperature 
and follows the peculiar scenario of Berezinskii, Koster- 
litz, Thouless [l|. On the contrary, the system un- 
dergoes Bose-Einstein condensation (BEC) only at zero 
temperature since long- wavelength phase fluctuations 
destroy long-range order. The analytical descriptions of 
2D systems include mean-field Gross-Pitaevskii (GPE) 
theory 0, beyond mean- field approaches @, Q, and 
numerical methods 



Recently, Bose-Einstein condensation was achieved in 
52 Cr atoms (lfjj . Chromium possesses a large perma- 
nent magnetic moment making the experimental obser- 
vation of effects of dipole-dipole interactions possible [11 1 . 
Moreover, by tuning its s-wave scattering length by a 
Feshbach resonance, one can realize a system with purely 
dipolar interactions. The first study of the Feshbach res- 
onances in 52 Cr atoms appeared recently (T^. As shown 
in Ref. [lj|, this procedure would permit to tune the ef- 
fective quasi-two-dimensional coupling constant g2D- In 
addition, the Bose-Einstein condensation of 2D dipolar 
excitons in quantum wells was recently observed in lu- 

In two coupled quantum 
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mmescence experiments 
wells, one containing only holes, and the other only elec- 
trons, holes and electrons might couple forming indirect 
excitons. Alternatively, the exciton dipole moment can 
be induced by normal electric field in a single quantum 
well. Spatial separation between hole and electron sup- 
presses recombination and greatly increases the lifetime 
of an exciton. If the separation between excitons are 
greater than the electron - hole separation D, an indirect 
exciton can be approximated as a boson with dipolar mo- 
ment oriented perpendicularly to the plane. 



A homogeneous system of N dipoles is described by 
the Hamiltonian 
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where M is the mass of a dipole. The coupling constant 
Cdd depends on the nature of the dipolar interaction. 
For cold atoms with large permanent magnetic moment 
m one has Cdd = Horn 2 , where hq is the permeability of 
free space. If the dipole moment is induced by an electric 
field E, the coupling constant is Cdd = E 2 a 2 /eo, where 
a is the static polarizability and eo the permittivity of 
free space. In the case of excitons Cdd = e 2 D 2 /e, where 
e is the electron charge and £ is the dielectric constant of 
the semiconductor. 

We suppose that the atoms are confined in a very tight 
pancake trap, so that dynamically the system is two- 
dimensional, and the dipole moments are aligned perpen- 
dicularly to the plane. This stabilizes the system since 
interactions are purely repulsive. 

The properties of Hamiltonian {T]) are governed by the 
two-dimensional gas parameter na 2 , n being the 2D den- 
sity and a the 2D s-wave scattering length. The limit 
na 2 — ► corresponds to a weakly interacting regime 
where the mean-field (MF) theory can be applied. The 
MF energy per particle is proportional to the coupling 
constant: E MF /N — gn/2. The leading contribution to 
the purely 2D coupling constant was first obtained b 
Shick[3 and rigorously derived for GPE by Lieb et al. 
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A peculiarity of a two-dimensional system is that the cou- 
pling constant depends on density, while it is independent 
of density in three- and one- dimensional systems. This 
factor causes particular energetic properties (for exam- 
ple, in the frequencies of collective oscillations). Calcula- 
tion of subleading corrections to MF ([2]) is a difficult and 
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long-standing problem and it was addressed by different 
authors [B[ giving sometimes controversial results. 

In a 2D system a Bose-Einstein condensate can only be 
formed at zero temperature. The system is completely 
condensed in the limit na 2 —>■ 0, while in a denser system 
the condensate density no gets depleted. The departure 
from the fully condensed state is given by the perturba- 
tive expression Q 
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For small momenta the excitation spectrum is linear 
Efc = h\k\c, with c being the speed of sound, which cor- 
responds to having phonons in the system. This defines 
a linear behavior of the static structure factor 
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in agreement with Feynman formula [l5j that relates 
the excitation spectrum to the static structure factor 
£ k = h 2 k 2 /{2MS k ). Excitations at very high momenta 
fc|/n 1/2 > 1 are free particles £ k = h 2 k 2 /2M and 
Sk = 1- The intermediate behavior of the static struc- 
ture factor at Ifcjn 1 / 2 1 can be approximated in the 
weakly-interacting regime by assuming the Bogoliubov 



excitation spectrum £ k — Jh k 2 pt/M + (Ti k 2 /2M) 2 , 



-iBog. 



y^Wnn/l In na 2 



< 1 



(5) 



We analyze the ground-state properties of the system 
described by the Hamiltonian (|TJ) resorting to the Diffu- 
sion Monte Carlo (DMC) technique. The guiding wave 
function is of Jastrow type, with a two-body correlation 
terms corresponding at short distances to a solution of 
the two-body scattering problem and smoothly matched 
to a long-range hydrodynamic asymptotics at larger dis- 
tances (further details in Ref. [7J). The parameters of 
the guiding wave function are optimized by energy min- 
imization in variational Monte Carlo calculations. The 
optimal guiding wave function is then used as an in- 
put to the DMC algorithm, which produces exact (apart 
from some statistical uncertainty) ground-state energies 
in a system of bosons. The pair distribution function 
and static structure factor are calculated using the tech- 
nique of pure estimators [l6| . The one-body density ma- 
trix is estimated using the extrapolated estimator us- 
ing DMC and VMC results [lf|. As finite-size effects 
are very important for the case of a dipolar interaction 
potential [7|, we perform calculations with system sizes 
N = 19; 25; 35; 50; 70; 100; 200 and do extrapolation of 
the energy to the thermodynamic limit N — > oo. 

DMC results for the ground-state energy as a function 
of the gas parameter na 2 are reported in Fig. Q] In the 
inset, we make a direct comparison to the mean- field en- 
ergy, Eq. We note that while the MF description is 




FIG. 1: Ground state energy per particle as a function of 
the gas parameter na 2 . Inset: in units of H 2 /(Ma 2 ), symbols 
- DMC data, line - mean- field energy Emf, Eq. ([2J. Main 
figure: in units of Emf, symbols - DMC data, solid line - 
best fit, Eq. ©, dashed line - mean- field energy. 



in a qualitative agreement with our data up to large val- 
ues of the gas parameter na 2 < 0.1, quantitatively there 
are noticeable differences. In order to test the quality of 
analytical equations of state we plot in Fig. Q]the energy 
in units of the MF energy. Even at the smallest consid- 
ered density na 2 = 10 -8 MF theory does not provide a 
quantitatively correct result. The behavior at such low 
densities is expected to be universal and defined by the 
gas parameter. Indeed, the energy of dipoles is almost 
indistinguishable from the energy of both soft and hard 



spheres [Uj when na < 10 . It is worth noticing that 



the situation is different from the 3D case where there 
is a perfect agreement with MF for values of the gas pa- 
rameter n^r>a^ D < 10~ 6 [l8(. In a ID system, the MF 
result is recovered with an accuracy better than 10~ 3 
already for nii)|ai£)| < 10~ 3 [3]. The reason for such 
a difference is that in three- and one- dimensional sys- 
tems the expansion comes in terms of the gas parameter, 
whereas in 2D the expansion is done in terms of 1/ In na 2 , 
and corrections to the energy include slowly converging 
terms In | In na 2 \ . This means that for experimentally 
interesting densities na 2 3> 10~ 8 , one must consider ad- 
ditional terms in the equation of state We tested 
a number of equations of state present in the literature 
[j| and none of them is able to reproduce correctly the 
DMC value of the energy even at na 2 — 10~ 8 . Thus, the 
only remaining way to obtain a quantitatively correct de- 
scription of a 2D dipolar gas is to perform a microscopic 
study. We fit the DMC equation of state in the density 
range 10 -6 < na 2 < 10 3 using the following function 




exp(a + ail + a2l 2 + a 3 l 3 + a i l i ), n > n c 
b\l + 62 In I + ^3 In In /) , n < n c 
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where / = |lnna 2 |. The best fit gives ao = 1.4552(4), 
ai = -1.3033(1), a 2 = 0.01267(2), a 3 = 7.7(2) ■ 10~ 5 , 
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FIG. 2: Condensate fraction as a function of the gas parame- 
ter na 2 , symbols - DMC data, dashed line - perturbative low 
density prediction, Eq. ©. 



FIG. 3: Static structure factor: symbols - DMC data, dashed 
lines - the low-momenta phononic behavior, Eq. Q, solid lines 
- Feynman prediction assuming Bogoliubov excitation spec- 
trum Eq. (0. Inset: pair distribution function. In descend- 
ing order of the slope at r = 0: na 2 = 9.4 • 10~ 9 ; 6.1 ■ 10~ 4 ; 
7.9 • 10" 2 ; 0.63; 10 



a 4 = -2.18(4) • lCT 5 , b = -2.225(20), h = 1, b 2 = 
3.294(5), 6 3 = -6.65(1), n c = 1.8 • 1(T 4 . Numbers in 
parentheses indicate the uncertainty of the fitting proce- 
dure. The resulting curve is shown in Fig. [1] 

The condensate fraction no/n is obtained from the 
non-diagonal long-range asymptotics of the one-body 
density matrix [7(. In Fig. [21 we plot results of no/n 
as a function of the gas parameter. We find that the 
predictions of Bogoliubov theory, Eq. ((3]), are in a good 
agreement in the weakly interacting regime. This should 
be contrasted with the energy, where the difference at the 
same density is significative. Nevertheless, at larger den- 
sities deviations are found. We note that at the densities 
of physical interest for excitonic systems na 2 « 1 — 10, 
the condensate is highly suppressed, making the applica- 
bility of the Gross-Pitaevskii equation, where the system 
is thought to be fully condensed, hardly possible. More- 
over, in this regime it is expected that the universal de- 
scription in terms of the s-wave scattering length breaks 
down and the properties of the system depend crucially 
on details of the potential. 

In Fig. [3j we plot the static structure factor Sk calcu- 
lated in a system with N ~ 200 particles in a wide range 
of the gas parameter 9.4 • 10 -9 < na 2 < 10. The static 
structure factor is obtained as a Fourier transformation 
of the pair distribution function. The low-momentum 
part of Sk is in agreement with the linear phonon pre- 
diction (j4|), shown in Fig. [3] with short-dashed lines. The 
slope is proportional to the speed of sound c, which, in 
Fig. [31 is calculated from the fit We see that for 

the two smallest considered densities Bogoliubov theory 
of weakly interacting Bose gas provides a relatively good 
description of the static structure factor. In the inset 
of Fig. [3] we show the characteristic behavior of the pair 
distribution function g2{f)- In the regime of low densi- 
ties <?2 (f) goes smoothly from zero (the dipole potential 



is highly repulsive when r — > 0) to the constant value 
.92(f) — > n 2 ,?! 1 / 2 )^ ^> 1, which is a typical behavior for 
a liquid (gas). We note that in this regime at small dis- 
tances the pair distribution function is proportional to 
the square of the two-body scattering solution /^(r): 
cj2(r) — const ■ |/ts(^)| 2 , while the coefficient of propor- 
tionality is a many-body property and can not be ex- 
tracted from the two-body problem. At large densities 
na 2 w 1 a peak at mean interparticle distance r = n~ 1 / 2 
starts to be formed, which shows that the system is more 
correlated in this regime. Indeed, it was found@,[i that 
the system of dipolcs experience a quantum phase tran- 
sition to a crystal if the density is greatly increased. 

Differences in the equation of state have consequences 
also in trapped systems, leading to different values of 
the release energy, the size of the cloud, etc. One of 
the most precise techniques, which can be used for test- 
ing the equation of state, is the measurement of the fre- 
quency of the "breathing" mode produced by a sudden 
change of the frequency of the harmonic trapping in the 
2D plane. In order to calculate the frequency of this 
mode in a trapped system we have used the local den- 
sity approximation (LDA) [2oT |. Within this approxima- 
tion all properties depend on a characteristic parameter 
A' r a 2 /a 2 , where ciho is the oscillator length of the confine- 
ment in the 2D plane. The frequencies of the "breathing" 
mode in a 2D system were calculated in Ref. [20( using 
scaling approach and assuming that the amplitude of os- 
cillations is small. In a symmetric trap this frequency is 
given by 

The mean square size of the cloud is calculated within 
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FIG. 4: Reduced square f2 2 ju\ of the lowest monopole fre- 
quency ("breathing" mode) in an isotropic trap, calculated 
as Eq. ©, for the DMC (solid line), MF (dash-dotted line) 
equations of states as a function of LDA parameter. For com- 
parison we plot n — > asymptotic Q = 2uJho (dashed line). 



local density approximation. This frequency provides a 
sensitive tool for testing the equation of state and can 
be measured in near future in experiments. An impor- 
tant conclusion of the present work is the failure of the 
commonly used Gross-Pitaevskii mean-field theory (dif- 
ferently from 3D and ID systems) to describe quanti- 
tatively energetic properties of a weakly interacting 2D 
system of dipoles. This failure is present already in the 
universal regime, i.e., where the shape of the interaction 
potential is no longer important, suggesting that special 
care should be taken in the description of properties of 
any 2D dilute Bose system. 

The work was partially supported by (Spain) Grant 
No. FIS2005-04181, Generalitat de Catalunya Grant No. 
2005SGR-00779 and RFBR. G.E.A. acknowledges post 
doctoral fellowship by MEC (Spain). 



LDA assuming that (r 2 ) 3> a| . This condition is easily 
satisfied if the number of dipoles in the trap N is large. 

In Fig. [U we show frequencies of the "breathing" mode 
f2 obtained for different equations of state. In the regime 
of weak interactions Na 2 /af — > the frequency goes 
asymptotically to a constant value fl = 2wh a . This value 
corresponds to the mean-field result /i = gn, with a cou- 



pling constant g independent of density [22j. Differently 
to three- and one- dimensional systems, the purely- two- 
dimensional regime of constant f2 can not be probably 
reached in experiments, since the system should be ex- 
tremely dilute. Indeed, even for the smallest considered 
values of the LDA parameter Na 2 /a\ — 1CP 5 the square 
of the frequency fi 2 » 4.2w 2 D is considerably different 
from the limiting one. Already for this value of the 
LDA parameter we find differences between the DMC 
and mean-field equations of state. This discrepancy be- 
comes more evident when the system enters into a regime 
of stronger interactions Na 2 ja 2 lo « 10 _1 . On the other 
hand, in the characteristic densities of excitonic experi- 
ments no 2 1 — 10 the mean-field description, Eq. 
fails completely and data should rely on DMC results. 
We note that differences of several percent can be re- 
solved in present high-precision experiments with cold 



gases (see, for example, Ref. [21 



In conclusion, we have performed a fully microscopic 
description of a weakly interacting system of dipoles us- 
ing the DMC method. The ground-state equation of 
state, as well as results for the static structure factor, 
pair distribution function and condensate fraction are re- 
ported. It has been shown that Bogoliubov theory pro- 
vides a good description of correlation functions in the 
dilute regime. For a harmonically trapped systems we 
calculate the frequency of the "breathing" mode using 
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